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Abstract 

In plants, a two component system (TCS) composed of sensor histidine kinases (HKs), histidine phospho- 
transfer proteins (HPs), and response regulators (RRs) has been employed in cyto!<inin signal transduction. 
ATCS exhibits important functions in diverse biological processes, includingplant growth, development, and 
response to environmental stimuli. Conducting an exhaustive search of the Chinese cabbage genome, a total 
of 20 HK(L) (1 1 HKs and 9 HKLs), 8 HP (7 authentic and 1 pseudo), and 57 RR (21 Type-A, 1 7 Type-B, 
4 Type-C, and 1 5 pseudo) proteins were identified. The structures, conserved domains, and phylogenetic 
relationships of these protein-coding genes were analysed in detail. The duplications, evolutionary patterns, 
and divergence of the TCS genes were investigated. The transcription levels of TCS genes in various tissues, 
organs, and developmental stages were further analysed to obtain information of the functions of these 
genes. Cytokinin-related binding elements were found in the putative promoter regions of Type-A BrRR 
genes. Furthermore, gene expression patterns to adverse environmental stresses (drought and high salinity) 
and exogenous phytohormones (tZand ABA) were investigated. Numerous stress-responsive candidate genes 
were obtained. Our systematic analyses provided insights into the characterization of the TCS genes in 
Chinese cabbage and basis for further functional studies of such genes. 
Key words: two-component system; Chinese cabbage; phylogeny; evolution; expression 



1 . Introduction 

Cytokinins are N ^-substituted adenine derivatives 
that have significant functions in various aspects of 
plant growth and development, including apical dom- 
inance, shoot or root branching, leaf expansion, lateral 
bud growth, photosynthesis, seed germination, floral 
transition, and leaf senescence.' Cytokinin signal trans- 
duction is mediated by a two-component system (TCS), 
which is similar to the two-component phosphorelay 
system that allows bacteria to sense and respond to 
environmental changes.^ Although this system was 



originally discovered in bacteria, two-component sig- 
nalling elements have also been identified in fungi, 
slime molds, and plants. However, a canonical histidyl- 
aspartyl phosphorelay is not found in animals.^"^ A 
simple TCS involves a histidine (His) sensor kinase and 
a response regulator (RR)^'^ His kinase (HK) perceives 
environmental stimuli via an input domain and autop- 
hosphorylates in a conserved His residue in a kinase 
domain. This phosphoryl group is then transferred 
to a conserved Asp residue on the receiver (Rec) 
domain of an RR. An RR undergoes phosphorylation, 
thereby modulating its ability to mediate downstream 
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signalling. Eukaryotes, such as yeast and plants, have 
evolved a more complex multi-step TCS with additional 
phosphorylation steps and His phosphotransfer (HP) 
proteins.^ Multi-component phosphorelay systems 
employ HK signal transduction in a multi-step His- 
Asp-His-Asp phosphotransfer process. 

Genes involved in the TCS have been extensively 
studied in Arabidopsis. The completion of the 
Arabidopsis genome sequence has revealed 55 genes 
encoding putative HK (AHK), HP (AHP), RR (ARR), and 
related proteins. Arabidopsis HKs are grouped into 
three mainsubfamilies:cytol<inin receptor; ethylene re- 
ceptor;and phytochrome. In addition, three HKs (CKI1 , 
AHK5, and AHKl) are present in Arabidopsis; such 
kinasesdo not belong in these groups. Theove rail struc- 
tures com prise a hybrid sensor HK containing a variable 
input domain, several N-terminal transmembrane 
domains, and a transmitter domain with a conserved 
structure, which includes the His residue that is the 
site of autophosphorylation, and a fused Rec domain.^ 
However, the transmitters of three ethylene receptors 
(ETR2, ERS2,and EIN4) and phytochromes lack key resi- 
dues in highly conserved HK sequence motifs; thus, 
these transmitters unlikely perform HK activity. For 
this reason, these molecules are called diverged HKs of 
two-co mponentelements.^°~^^ Mo reove r, t h ree cy to- 
kinin receptors (AHK2, AHK3, and AHK4) share a 
cyclase/HK-associated sensory extracellular (CHASE) 
domain, which is a putative recognition site for cytoki- 
nins.^'^ Five ethylene receptors (ETRl, ERS1, ETR2, 
ERS2, and EIN4) also share an ethylene-binding 
domain (C2H4). Five phytochromes (PHYA, B, C, D, 
and E) are characterized by a chromophore-binding 
domain (PHY) and two PAS (Per/Arndt/Sim) folds. 
AHPs, which contain a highly conserved xHQxKGSSxS 
motif, can mediate the transfer of a phosphate group 
from the Rec domain of an AHK to the Rec domain of 
an ARR.^^ AHP6 is a pseudo-His-containing phospho- 
transfer (HPt) protein without a conserved histidine 
residue; thus,AHP6 cannot function asa phosphotrans- 
fer protein. In fact,AHP6 functions as a negative regula- 
tor of cytokinin signalling by inhibiting phosphorelay 
from phosphorylated AHP1 to ARRI.^^ The canonical 
members of ARRs can be divided into two major 
classes: Type-A ARRs and Type-B ARRs, according to 
their amino acid sequences and conserved domains. 
Type-A ARRs are primary cytokinin response proteins 
comprising a Rec domain along with short C-terminal 
extensions.Type-BARRsarecomposedof an N-terminal 
Rec domain and a large C-terminal output domain with 
a Type-B signature GARP (GOLDEN/ARR/Psrl ) motif of 
~ 60 amino acids ' ^ that are distantly related totheMyb 
DNA-binding superfamily.^ ^'^ ^ Type-C ARRs are also 
present, and the members of this protein class were ori- 
ginally considered as Type-A ARRs because these mole- 
cules have the same domain structure.^" In contrast to 



the Type-A ARR genes, Type-C ARR genes cannot be 
induced by cytokinin,^°~^^ and their functions in the 
cytokinin signalling pathway remain unknown. In 
addition, diverge RRs, referred to as pseudo-RRs (PRRs), 
lack the conserved Asp for phosphorylation.^^ As a 
highly characterized pseudo-RR, Clock PRR contains a 
distinct plant-specific Co, Col and Tocl (CCT) motif in 
their C-terminal extensions and functions in the regula- 
tion of circadian rhythms.^"^"^^ 

Abiotic stress conditions, such as drought and high 
salinity, are two of the most common stress factors 
that adversely affect plant growth and yield. Plants 
have evolved a complex signalling network at molecu- 
lar, cellular, and system levels to survive and flourish 
in varied environments.^^ Phosphorylation, which is 
mediated by TCS genes, is a key mechanism of stress 
signal transduction in plant cells.^^ Increasing evi- 
dence demonstrates that TCS genes are involved in 
the response pathway to environmental stimuli 
in Arabidopsis, rice, and soybean. ^^"^^ To date, TCS 
genes have been identified at a genome-wide scale 
in various plant species, including Arabidopsis,^ ^'^^ 
j.j^g^26,38-4i Lotus japonicus,'^^ soybean, ^^'^^ 
maize,"^^'"^"^ and Physcomitrella patens,'^^~'^^ among 
others. However, the knowledge of their potential 
functions in stress adaptations remains confined to 
Arabidopsis and rice. Studies on TCS genes in Chinese 
cabbage are yet to be conducted. 

Chinese cabbage {Brassica rapa ssp. pel<inensis) is one 
of the most important B. rapa crops and considered as 
an economically important vegetable worldwide 
because of its high yield and good quality. The whole- 
genome sequencing of 6. rapa (Chiifu-401 -42) by 
the Brassica rapa Genome Sequencing Project 
Consortium'^^ enables us to undertake a genome- 
wide identification and functional analysis of the gene 
families related to the morphological diversity and 
agronomic traits of Brassica crops. Furthermore, the 
'A' genome of 6. rapa is an important resource to 
study the evolution of polyploidy genomes and poten- 
tial strategies to improve 6rass/cc/-related crops genetic- 
ally.^° Given the significance of TCS genes in diverse 
biological and physiological processes, including our 
main interest area of abiotic stress responses, whole 
genome-wide analysis was performed to identify 
the HK, HP, and RR proteins involved in TCS in Chinese 
cabbage. The gene structures, conserved domains, 
and phylogenetic relationships of the TCS genes were 
analysed in detail. A comprehensive analysis of tissue 
or organ-specific expressions of TCS genes in Chinese 
cabbage, and their expression profiles to abiotic stresses 
(drought and high salinity) and exogenous phytohor- 
mones (tZ and ABA) were conducted. Moreover, the 
gene duplications and evolutionary patterns of 
the TCS genes in Chinese cabbage and Arabidopsis 
were determined. Our results provided a framework of 
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future functional analyses of the TCS genes and lay a 
good foundation for the utilization of potential genes 
used for breeding to enhance plant production, 
quality, and stress-resistance in Chinese cabbage. 



2. Materials and methods 

2.7. Identification of TCS genes in Chinese cabbage 

Arabidopsis TCS protein sequences were used as seed 
sequences to search the Brassica Database (BRAD) 
Version 1.1 (http://brassicadb.org/brad/)'^^'^^'^^ and 
the NCBI database (www.ncbi.nlm.nih.gov). The search 
was based on a BLASTP search with an expected value 
of 1 00. The protein sequences of the identified TCS 
members were used asqueriesto reconfirm the multiple 
databases to ensure that no additional related genes 
were missing from the database. All of the sequences 
that satisfied the requirements were analysed by using 
the Pfam database (http://pfam.janelia.org/), the 
SMART database (http://smart.embl- heidelberg.de/), 
and the Conserved Domain Database of the NCBI 
(http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb. 
cgi)^^ to eliminate the genes that did not contain the 
known conserved domains and motifs of the TCS 
members. The genomic information concerning on 
the chromosome locations of the TCS genes and the 
amino acid sequences of TCS proteins were obtained 
from BRAD. HKL members in soybean were identified 
and supplied as described. ^^Sequence identities were 
calculated using the MatGATsoftware.^^ 



2.2. Gene structure, motif recognition, multiple- 
sequence alignment, and phylogenetic analyses 
The gene structure schematic of the TCS genes 
was illustrated using the Gene Structure Display 
Server (http://gsds.cbi.pku.edu.cn/). To identify the 
conserved motifs of the TCS protein sequences, we 
used the online MEME (http://meme.sdsc.edu/ 
meme/meme.html).^^ The deduced amino acid 
sequences of the conserved HK domain, Rec domain, 
Hpt domain, Myb domain, and CCT motif were 
aligned usingtheClustalXprogram.^^ Forthe phylogen- 
etic relationships of HK, HP, and RR genes in B. rapa 
and Arabidopsis thaliana, the identified HK, HP, and RR 
protein sequences were separately aligned with a 
gap opening penalty of 10 and a gap extension 
penalty of 0.2 by using ClustalW implemented in the 
MEGA5.0 software (http://www.megasoftware.net/). 
The phylogenetic trees were further constructed using 
the neighbour-joining (NJ) method, in which Poisson 
correction, pairwise deletion, and bootstrapping 
(1 000 replicates; random seeds) were considered as 
parameters. 



2.3. Analysis of the putative promoter regions of the TCS 
genes in Chinese cabbage 

The upstream sequences (1 000 bp) of the transcrip- 
tional start site of each TCS gene were chosen as the 
putative promoter regions to identify the abiotic stress- 
related and phytohormone-related c/s-elements. The 
PlantCARE website (http://bioinformatics.psb.ugent. 
be/webtools/plantcare/html/search_CARE.html) was 
used to identify the putative c/s-regulatory elements 
along the promoter sequences.^" Besides, the recently 
reported stress-responsive c/s-motifs that were not 
covered,^^-*^^ such as ICErl, ICEr2, NACRS, ZFHDR, 
MYBR, MYCR, CRT, ABRE, EE, G-box, T/G box, and CE3, 
werealsoused asqueriestosearchagainstthe promoter 
sequences. In addition, the upstream sequences 
(1 000 bp) of the Type-A RR genes were used to query 
the PLACE database (http://www.dna.affrc.go.jp/ 
PLACE/)^^ and search for the three regulatory motifs: 
the cytokinin-responsive Type-B ARRl -binding ele- 
ments AGATT*^*^ and GATCTT' and the cytokinin- 
enhanced protein-binding element TATTAG.^° 

2.4. Chromosomal localization, gene duplications, and 
evolutionary analysis of the TCS genes in Chinese 
cabbage and Arabidopsis 

The Brass/cfl Genome Browse (http://brassicadb.org/ 
cgi-bin/gbrowse/cbgdbl 1 /) was used to map the posi- 
tions of the TCS genes in the physical maps of 1 0 6. rapa 
chromosomes. The TCS genes of A. thaliana were dis- 
tributed and visualized using a chromosome map tool 
(http://www.arabidopsis.org/jsp/ ChromosomeMap/ 
tool.jsp). Tandem duplications were defined whether 
or not the two genes are separated by four or fewer 
gene loci.^^ Segmental duplications were identified by 
synteny analysis using an online tool (http://chibba. 
agtec.uga.edu/duplication/). ''^ The synteny analysis 
of the TCS genes was performed using PGDD (http:// 
chibba.agtec.uga.edu/duplication/). The occurrence 
of duplication events and homologous genes diver- 
gence as well as selective pressure on duplicated genes 
were estimated by calculating synonymous (Ks) and 
non-synonymous substitutions (Ka) per site between 
the duplicated gene pairs by using the CodemI proced- 
ure of the PAML program. The divergence time was 
calculated at a neutral substitution rate of 1 .5 x 1 0~^ 
substitutions per site per yearfor the chalcone synthase 
gene (Chs).'''^ 

2.5. Plant growth, treatments, and sampling 
Brassica rapa subsp. pekinensis line Chiifu-401 -42 

was grown in an experimental farm in Zhejiang 
University. The roots, floral stems, leaves, little buds 
(<1.6mm), medium-sized buds (1.6-2.8 mm), big 
buds (>2.8 mm), flowers, sepals, petals, stamens, 
pistils, little siliques, medium-sized siliques, and big 
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siliques of the plant were sampled to analyse the tissue- 
or organ-specific expressions. Methods and sites of 
sampling were described previously.^^ Besides, little 
siliques, medium-sized siliques, and big siliques were 
defined as siliques at 3, 1 4, and 2 5 days after pollination 
and fertilization. 

All of the seedlings used for treatment were grown at 
25 ± rc for ~3 weeks with a 1 6-h light/8-h dark 
photoperiod. The nutrient solution was supplied with 
200 mM NaCi for salt treatment. The total roots were 
separately collected at 0, 3, 8, 24, and 48 h after 
stress induction. Three-week-old seedlings were with- 
held from watering to simulate drought conditions, 
and the seedlings were divided into four levels (0, 1, II, 
and III) based on the degree of the drought symptoms 
as previously described. The total roots were separ- 
ately collected at these four levels. For mock control, 
seedlings were grown at normal conditions with nutri- 
ent solution and samples were collected in the mean- 
while. For plant hormone treatment, 3-week-old 
seedlings were sprayed with the following: 1 00 jxM 
trans-zeatin (tZ) for the cytokinin treatment and 
100 |xM abscisic acid (ABA) for the ABA treatment. 
Only the second true leaves were sampled to minimize 
differences. The leaves were sampled at 0, 0.5, 1 , 2, and 
4 h after spraying, and the control sample was sprayed 
with double distilled water alone. All the samples were 
immediatelyfrozen in liquid nitrogen and stored in a re- 
frigerator at -7 5°C. 

2.6. RNA extraction and quantitative reverse 
transcription PCR analysis 
RNA extraction, reverse-transcription, and quantita- 
tive reverse transcription PCR (qRT-PCR) reactions were 
performed as previously described^^ using the primers 
listed in Supplementary Table SI . The specificity of the 
reactions was verified by melting curve analysis, and 
the products were furtherconfirmed by agarose gel elec- 
trophoresis. The BrCyp gene was used as the reference 
gene to study different organs or tissues and develop- 
mental stages; the BrUBC30 gene was used in abiotic 
stress and plant hormone treatment studies as recom- 
mended by Xiao etal/^The comparative AAQ method 
was used to calculate the relative expression levels of 
different genes. qRT-PCR results were clustered using 
the average linkage method with Pearson correlation 
distance metric by Multiple Array Viewer.^'' 

3. Results and discussion 

3.1 . Identification and annotation of the TCS genes 
in Chinese cabbage 
A total of 85 TCS gene members were identified 
based on BLASTP search results against the BRAD and 
NCBI databases from the 6. rapa genome. The DNA, 



CDS, protein sequences, and promoter regions of all 
the identified TCS members were provided in 
Supplementary Dataset 1. A total of 2 0 6rHK(L) pro- 
teins were divided into two groups, BrHKs and BrHKLs, 
according to the characteristics of their conserved 
domains (Supplementary Table S2). Eleven BrHKs con- 
tained a conserved histidine-kinase domain (HK), 
whereas nine BrHKLs comprised a diverged histidine- 
kinase-like domain (HKL). Seven BrHPs with conserved 
HPt domain and BrPHPI with a pseudo-HPt lacking 
the His phosphorylation site were retrieved 
(Supplementary Table S3). A total of 57 RR proteins 
were further classified into 21 Type-A RRs, 1 7 Type-B 
RRs, 4 Type-C RRs, and 1 5 pseudo-RRs based on their 
conserved domains and motifs (Supplementary Table 
S4). 6rRR2 2 lacked the Myb domain, but 6rRR2 2 was 
still assigned to Type-B RRs because this gene exhibited 
high similarities with the proteins in this group. The full- 
length gene was unlikelyobtained;thus,furtheranalysis 
should be conducted to obtain the full-length clone and 
to classify this gene in the Type-B RR subfamily. Many 
TCS proteins in Chinese cabbage shared high sequence 
identities with thek Arabidopsis counterparts, indicat- 
ing the conservative evolution of TCS genes. To show a 
clear and systematic understanding of the number of 
TCS genes in plant genomes, we summarized the TCS 
genes that have already been identified in several 
plants (Table 1). Considering that Brassica genomes 
have undergone another whole-genome triplication 
(WGT) after speciation from A. thaliana, we presumed 
that the TCS genes in Chinese cabbage should have ap- 
proximately three times as many members as that of 
Arabidopsis. However, the TCS genes identified in 
Chinese cabbage were only 1.55 times as many as 
those found in Arabidopsis. This result suggested that a 
substantial loss of genes after a hexaploid was formed 
by WGT. This finding is consistent with that in previous 
studies, in which a total of 41 174 protein-coding 
genes have been identified in the 6. rapa genome; 
such genes were roughly 1 .5 times as many as those 
found in A. thaliana (27 A 1 1 genes in TAIRIO).'^'' A 
total of 21 and 18 HK genes are found in the 
genomes of Glycine max and P. patens; however, this 
finding has raised a question on the numerous HK 
genes present and the mechanism by which these HK 
genes were obtained during evolution. The origin and 
evolutionary pattern of the TCS signalling pathway 
were investigated by surveying the genomes of several 
sequenced key plant species ranging from unicellular 
algae, moss, and lycophytes to higher land plants, in- 
cluding/4rafa/c/ops/s and rice, to determine the proteins 
involved in cytokinin signal transduction.'^^ CHASE 
domain-containing proteins were not detected in the 
genomes of any of the algal species. Physcomitrella 
patens was possibly one of the most basal plant 
species containing all of the components necessary 
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Species 


HK 


HPt 


Type-A RR 


Type-B RR 


Type-C RR 


Pseudo-RR 


Total 


References 


Arabidopsis thaliana 


8 


6" 


10 


12 


2 


9 


47 


1 5 


Oryza satiua 


8 


5" 


13 


13 


2 


8 


49 


41 


Lotus japonicus 


14 


7 


7 


1 1 


1 


5^ 


40 


42 


Glycine max 


21 


13 


18 


1 5 


3 


13 


83 


28 


Zea mays 


1 1 


9'= 


16 


9 


3 




59 


43,44 


Physcomitrella patens 


1 8 


3 


7 


5 


2 




39 


45-47 


Brassica rapa 


1 1 


8' 


21 


1 7 


4 


1 5 


76 


This work 



^Five authentic and one pseudo-HPts. 
''Two authentic and three pseudo-HPts. 
'^Seven authentic and two pseudo-HPts. 
'^Seven authentic and one pseudo-HPts. 
'^Only clock-associated. 



for cytokinin signalling. Moreover, the number of the 
cytokinin receptor genes did not increase in the same 
manner as that of the other components. This result is 
also true for the cytokinin receptor genes of Chinese 
cabbage in ourstudy. In contrast,Chu etal.'^^ identified 
1 1 ZmHKs with highly conserved CHASEdomains. Eight 
ZmHKs comprised three conserved domains (CHASE, 
transmitter, and Rec domains). Three ZmHKLs con- 
tained a CHASEdomain, but these proteinsdo not com- 
prise an integral transmitter domain or Rec domain, 
suggesting that the CHASE domain is highly conserved 
in the evolution of ZmHKs. This surprising anomaly 
increased the complexity of the evolution of CHASE 
domain-containing cytokinin receptor genes in 
plants. By comparing the numbers of Type-A RR and 
Type-B RR genes in various species, it is noted that 
there was not always more Type-A RR genes than 
Type-B RR genes. Until now, we only have some knowl- 
edge that Type-B RR genes were already found in the 
unicellular algae while the Type-A RR genes, as the 
youngest subgroup of RR genes, only first appeared in 
the land plant species.'^^ Since RR members were iden- 
tified in only several species, more rules might be found 
if more RR genes were identified. 

3.2. Gene structure, conserved domain, 
and phylogenetic analysis 
To analyse the structural characteristics and con- 
served regions of the TCS genes, we mapped the gene 
structures containing exons and introns. We also exam- 
ined their conserved regions and motifs; furthermore, 
their putative protein sequences were aligned 
(Supplementary Fig. SI and Table S4). The HK, HPt, 
Rec, and Myb domains as well as the CCT motif of 
the two-component elements were well conserved. 
However, a group of diverged two-component ele- 
ments, such as BrHKLs, SrPHPl, and pseudo-SrRRs, 
lacked the conserved phosphorylation sites. In addition. 



to identify subgroups and reveal the evolutionary rela- 
tionships of the TCS genes in Chinese cabbage, 
Arabidopsis, rice, and soybean, the complete amino 
acid sequences of HK(L), HP, and RR proteins were used 
to perform multiple alignments and construct phylo- 
genetic trees. Besides, 1 5 GmHKLs were identified and 
supplied for a more complete phylogenetic tree 
(Supplementary Table S5 and Dataset 2). Figure 1 
shows that HK(L) members were further divided into 
six subfamilies: cytokinin receptor subfamily; AHK5/ 
CKI2 subfamily;AHK1 subfamily;CKIl subfamily; ethyl- 
ene receptor subfamily; and phytochrome subfamily. 
Cytokinin receptor subfamily comprised four BrHK 
proteins (6rHK6-9) with HK, HATPase, Rec, transmem- 
brane (TM), and CHASE domains. The CHASE domain is 
specific for cytokinin receptor genes and necessary for 
cytokinin binding; this result indicated that Brl-ll<6-9 
genes may function as cytokinin receptors. The pres- 
ence of a TM domain in the cytokinin receptor genes 
of Chinese cabbage indicatestheendomembrane loca- 
tion of cytokinin-binding sites.^*''^'' AHK5/CKI2 sub- 
family comprised three SrHK proteins (6rHK2-4) with 
HK as well as HATPase and Rec domains; however, this 
subfamily lacked a TM domain. Both AHKl and CKIl 
subfamilies contained only one BrHK member with 
HK, HATPase, Rec, and TM domains. Ethylene receptor 
subfamily could be further divided into ethylene recep- 
tor subfamily I and ethylene receptor subfamily II. 
Ethylene receptor subfamily I was composed of two 
BrHK proteins (firHKI 0 and BrHKl 1) with TM, C2H4, 
GAF, and HK domains. SrHKl 1 lacked the Rec domain. 
Ethylene receptor subfamily II was composed of four 
BrHKL proteins (BrHKLI -4) with TM, C2H4, GAF, and 
HKL domains. Both SrHKLS and 6rHKL4 lacked the Rec 
domain. 6rHKL4 contained an additional MATH 
(meprin and TRAP homology) domain, which belonged 
to the TRAF-like superfamily. Based on the gene struc- 
ture and conserved domains analyses, we proposed 
that BrHI<L4 was formed by the fusion of two genes 
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Figure 1. Phylogenetic relationships of HK proteins and related proteins in Chinese cabbage, Arafo/dops/s, rice, and soybean. The bar indicates 
the relative divergence of the sequences examined. This figure appears in colour in the online version DNA Research. 
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(Supplementary Fig. S5A). However, the actual fusion 
mechanism remains unknown. C2H4 is an ethylene- 
binding domain, suggesting its potential function in 
ethylene receptors. Phytochrome subfamily contained 
five firHKL proteins (BrHKL5-9) with HKL, HATPase, 
GAP, PHY, and PAS domains. These BrHKLs lacked the 
Rec domain. PHY domain is a chromophore-binding 
domain. The presence of chromophore-binding 
domains demonstrated that they were candidate 
genes involved in phytochromes. The phylogentic tree 
constructed with all the authentic and pseudo-HPts of 
Chinese cabbage, Arabidopsis, rice, and soybean indi- 
cates the presence of subgroups (Fig. 2). BrHPl -4 and 
6rHP7 showed a close relationship with AHPl, AHP2, 
AHP3, and AHP5, which act as positive regulators in CK 
signalling.^" 6rHP5 and 6rHP6 were clustered into a 
same groupwithAHP4, which showed relativelydistinct 
genetic relationships from the other AHPs and might 
play negative roles in CK signalling.^" Besides, SrPHPI 
was close to AHP6, a pseudo-HPt that lacked the His 
phosphorylation site, acting inhibitory roles in CK 
responses by competing with other AHPs.^ ^ According 
to the classification method proposed by Schaller and 
co-workers,^ ^'^^ RR proteins were classified as Type-A, 
Type-B, Type-C, and pseudo-RRs (Fig. 3). All the Type-A 
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Figure 2. Phylogenetic relationships of HPt proteins and related 
proteins in Chinese cabbage, Arabidopsis, rice, and soybean. The 
bar indicates the relative divergence of the sequences examined. 



Figure 3. Phylogenetic relationships of RR proteins and related 
proteins in Chinese cabbage, Arabidopsis, nee, and soybean. The 
bar indicates the relative divergence of the sequences examined. 
This figure appears in colour in the online version of DNA Research. 
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RRs could be well clustered into a group, which com- 
prised 21 Chinese cabbage members (6rRR1-21) 
with one conserved Rec domain (Rec). Type-B RRs 
seemed to be more diverse with five subgroups (l-V). 
Seventeen Chinese cabbage Type-B RR members 
located in Type-B I RRs (erRR22-32), Type-B II RRs 
(6rRR36-38), and Type-B III RRs (erRR33-35). 
Notably, several OsRRs alone occupied two extra sub- 
groups. OsRR28,OsRR29,OsPRRl l,and OsPRRI 2 com- 
posed a single subgroup, Type-B IV and also OsRR31 -33 
composed another single subgroup, Type-B V. We sup- 
posed that this phenomenon might be due to the fact 
that rice, as a monocotyledon, had distant relationships 
with dicotyledons. Type-C RRs contained four Chinese 
cabbage members (6rRR39-42) with conserved Rec 
domain (Rec) alone. Pseudo-RR members were mainly 
divided into two groups, i.e. Clock PRR and Type-B PRR. 
These pseudo-RR members had relatively closer phylo- 
genetic relationships with Type-B RR members. Ten 
ClockBrPRR members fBrPRRI -1 O)containedthecon- 
served CCT motif and conserved pseudo-Rec domain, 
which lacked the conserved amino acid that is essential 
for phosphorylation. Five Type-B SrPRR members were 
located in two clades, in which one clade contained 
four Type-B firPRR members (BrPRRI 1-14) with the 
conserved pseudo-Rec domain and Myb domain; 
whereas the other clade contained BrPRRI 5 with the 
pseudo-Rec domain alone. The gene fusion phenom- 
enon was also observed. Besides, the pseudo-Rec and 



Myb domains, BrPRRI 0 also had a CDF (cation diffusion 
facilitator) domain. Based on the gene structure and 
conserved domain analyses, BrPR/? 7 0 was also proposed 
to be formed byfusion of two genes (Supplementary Fig. 
S5 B). They are not alone, in fact, we previously identified 
that BrC/0C7 -3 was formed by the fusion of the C/OCand 
TPR gene.^^ 

On the whole, all the TCS members from Chinese 
cabbage, Arabidopsis, rice, and soybean were classified 
and grouped well in our phylogenetic trees, which was 
consistent with previous results.^^'^^ However, it is 
evident that TCS members from the same species were 
interspersed in each subgroup, especially the mono- 
cotyledonous rice, suggesting that the closer phylogen- 
etic relationships of the species, the closer phylogenetic 
relationships of homologous genes. For Chinese 
cabbage and Arabidopsis, each gene of Arabidopsis 
almost corresponded to 1-3 homologous Chinese 
cabbage genes in each clade of the phylogenetic tree, 
further supporting that Chinese cabbage genome 
underwent another WGT after divergence from 
Arabidopsis and this event occurred after Brassicaceae 
diverged from Gramineae and Leguminosae. 

3.3. Chromosomal distribution and duplications of TCS 
genes in Chinese cabbage and Arabidopsis 
Chromosomal location of each TCS gene was deter- 
mined from the genomic sequences of Chinese 
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Figure 4. Chromosomal mapping of TCS genes in Chinese cabbage. The arrows next to gene names show the transcription direction. 
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cabbage; in addition, 83 TCS genes were mapped on the 
1 0 chromosomes with an apparently uneven distribu- 
tion (Fig. 4). The two unmapped genes, BrHKG and 
BrRR40, were located on ScaffoldOOOl 04 and 
ScaffoldOOOl 91, respectively. Gene duplication events 
are important to gene family evolution because dupli- 
cated genes provide the raw materials for the generation 
of new genes, which in turn facilitate the generation of 
newfunctions.^^ Brassicaceae genomes have undergone 
three rounds of whole-genome duplication; these 
genomes are referred to as 1 R, 2R, and 3R, which are 
equivalent to the -y, |3, and a duplication events; further- 
more, Brassica genomes have undergone another WGT 
(4R) after speciation from A. thaliana,^^~^'^ which led 
to significantly increased duplicated gene numbers in 
6. rapa. In plants, gene numbers are mainly expanded 
by segmental and tandem duplications in gene fam- 
ilies.^^ We examined tandem duplications and segmen- 
tal duplications of TCS genes in Chinese cabbage. As 
shown in Table 2, multiple TCS genes were involved in 
segmental duplications, whereas only one tandem 
duplicated gene pair was identified; therefore, the ex- 
pansion of TCS genes in Chinese cabbage was mainly 
attributed to segmental duplications. Compared with 
duplicated TCS genes in Chinese cabbage, TCS genes in 
A. thaliana were also involved in segmental duplications, 
but with smaller percentage. However, all the 1 OType-A 
ARR genes and half of the Type-B ARR genes were striking- 
ly segmentally duplicated. 

The synonymous substitution rates (Ks) and non-syn- 
onymous substitution rates (Ka) are measured to 
explore the gene divergence mechanism after duplica- 
tion. Large-scale duplication events are defined as sim- 
ultaneous gene duplication. Assuming a molecular 
clock, the Ks of these duplicates are expected to be 
similar over time. However, substantial rate variations 
are noted among genes.^^ To elucidate the evolution- 
ary patterns and divergence of the segmental dupli- 
cated TCS genes in Chinese cabbage and Arabidopsis, 
the/<sand/<c/ modes forthesegmentalduplicated para- 
logs and orthologs were determined (Supplementary 
Table S6). We used the relative Ks measure as the substi- 
tute for time to evaluate the divergence time between 
Chinese cabbage and Arabidopsis. The frequency 



distributions of the relative /<s values that were obtained 
from duplicated orthologous gene pairs, between the 
Chinese cabbage and Arabidopsis genomes and dupli- 
cated paralogous gene pairs in the Chinese cabbage 
and /4ral7/(yo/js/s genomes, are shown in Fig. 5. The rela- 
tive Ks distribution of the duplicated paralogous gene 
pairs for TCS genes in A. thaliana peaked from 0.8 to 
0.9, which indicated that these duplicated genes oc- 
curred at ~53 to 60 Mya, which corresponds to the 
3R time.^^'^^ The relative Ks distribution of the dupli- 
cated paralogous gene pairs for TCS genes in Chinese 
cabbage displayed peaks from 0.3 to 0.4 and from 0.8 
to 0.9, which suggests an additional duplication event 
occu rred at ~ 2 0 to 2 6 Mya after 3 R, wh ich corresponds 
to the 4R time.^^ Our results further confirmed the 
hypothesis that Brassica genomes have undergone 
another WGT after speciation from /4. thaliana. The ex- 
istence of gene duplications in Chinese cabbage and 
Arabidopsis that occurred during 3 R also demonstrated 
their conservation during the long-term evolution. For 
the duplicated orthologous gene pairs between 
Chinese cabbage and Arabidopsis, the relative /<s distri- 
bution showed a major peakfrom 0.4 to 0.5, which sug- 
gests that the two genomes were separated at ~2 7 to 
33 Mya; this finding was in agreement with previous 
studies that Brassica and Arabidopsis genomes diverged 
after the 3R event.^^"^^ In addition, another less 
obvious peak appeared from 0.7 to 0.8, which was re- 
sponsible for the 3R event before the speciation 
between Brassica and Arabidopsis genomes. 

3.4. Analysis of the putative promoter regions of TCS 
genes in Chinese cabbage 
c/s-Regulatoryelements, which are located upstream 
of genes and act as binding sites for TFs, have essential 
functions in determining the tissue-specif icor stress-re- 
sponsive expression patterns of genes.^^ To further 
understand transcriptional regulation and potential 
functions of these TCS genes, 1000 bp regions up- 
stream of the transcriptional start site were applied to 
identify c/s-regulatory elements. A number of abiotic 
stress-related (e.g. drought, high salinity, extreme tem- 
peratures, and wound) and hormone-related (e.g. ABA, 



Table 2. Comparison of duplicated TCS genes in Chinese cabbage and Arabidopsis 



Gene family 


Number of genes involved 
in segmental duplications 


Percentage (%) 


Number of genes involved 
in tandem duplications 


Percentage (%) 


AHK(L)/BrHI<(L) 


2/8 


1 2.5/40 


0/0 


0/0 


AHP/BrHP 


2/5 


33.3/62.5 


0/0 


0/0 


Type-A ARR/ BrRR 


10/21 


100/100 


0/0 


0/0 


Type-B ARR/ BrRR 


6/1 1 


50/64.7 


0/2 


0/1 1.8 


Type-CARR/BrRR 


0/4 


0/100 


0/0 


0/0 


Pseudo-ARR/BrRR 


0/12 


0/80 


0/0 


0/0 
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auxin, ethylene, GA, MeJA, and SA) c/s-elements were 
found in the putative promoters of TCS genes in 
Chinese cabbage. These c/s-elements are counted and 
classified in Supplementary Table S7. The occurrences 
of these c/s-elements suggested that these TCS genes 
might have potential functions in abiotic stress adapta- 
tions and various hormone signalling. 

GARP is a Type-B signature; and the binding sites of 
GARP are present in the Type-A ARR gene promoters, 
which suggest that Type-B ARR proteins may bind to 
these promoters and activate transcription.^^ Previous 
research also confirmed that cytokinin-dependent in- 
duction of the Type-A RR genes is partially dependent 
on transcriptional regulation by Type-B RR genes.^^'^^ 
We investigated the cytokinin-related bindingelements 
in the Type-ABr/?/? promoter regions. Cytokinin-respon- 
sive Type-B ARRl -bindingelements AGATTand GATCTT 
as well ascytokinin-enhanced protein-binding element 
TATTAG were widely evident (Supplementary Table S8); 
this finding indicated that Type-A BrRR genes might be 
transcriptionally regulated by Type-B BrRR proteins 
and cytokinin. Two ARRl -binding elements were 
found in some Type-A BrRR gene promoters; thus, we 
speculated that the TCS is conserved between Chinese 
cabbage and Arabidopsis at the level of system regula- 
tion, as well as gene structure. 

3.5. Expression profiles of TCS genes in various tissues, 
organs, and developmental stages 
Tissue- or organ-specific and developmental stage- 
related expression data are useful in identifying genes 



that are involved in defining the precise nature of indi- 
vidual tissues in a given developmental stage.^^ A clear 
expression profileofgenesof interest is an indispensable 
step to find and utilize agriculturally important genes. To 
obtainfirst insight into thef unctions of TCS genesduring 
the Chinese cabbage vegetative and reproductive devel- 
opment, qRT-PCR was used to analyse the transcription 
levels of TCS genes in varioustissues, organs, and devel- 
opmental stages. Totally, only 31 of the 61 TCS genes 
that could be distinguished by qRT-PCR were selected 
for expression analysis because of the high similarities 
of the TCS sequences, and the remaining 9 BrHKL 
genes and 1 5 SrPRR genes were not within the scope 
of our study. Based on Fig. 6 and Supplementary Fig. 
S6, data analyses indicated high variability in transcript 
abundance of the TCS genes in Chinese cabbage. A ma- 
jority of TCS genes showed relatively high expression 
levels in roots, such as BrHK3, BrHK4, BrHP1 , BrPHPl , 
BrRR3 , BrRR 12, and BrRR3 7 . Th is f ind i ng m ight be corre- 
lated with the fact that cytokinins are mainly synthe- 
sized in roots. Notably, BrRRIZ, BrRRZ, BrRR3, BrRR12, 
BrHK9, BrHPZ, and BrRR3 1 exhibited similar expression 
patterns and were grouped together. These genes were 
expressed abundantly in vegetative organs, floral buds, 
and developing siliques; however, these genes showed 
low transcripts in sepals, petals, and stamens; these find- 
ings indicated that they have spatio-temporal expres- 
sion characters. BrRRI 9 had a relatively higher 
expression level in floral buds and four whirl flower 
organs; this finding suggests its function in flowerdevel- 
opment. Genes with high transcripts in developing 
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Figure 6. Hierarchical clustering and heat map representation for the 
tissue- or organ-specific TCS gene expression profiles in Chinese 
cabbage. R, roots, S, floral stems, L, leaves, LB, little buds, MB, 
medium-buds, BB, big buds, F, flowers, Se, sepals, Pe, petals, St, 
stamens. Pi, pistils, LS, little siliques, MS, medium-sized siliques, 
and BS, big siliques. The expression levels of genes are presented 
using fold-change values transformed to Log2 format. This figure 
appears in colour in the online version of DMA Research. 



siliques may have essential functions in seed qualityand 
yield. We found that BrRR33 and BrRR38 were almost 
expressed exclusively in developing siliques, which 
demonstrates that they were candidate genes in im- 
proving seed quality and yield. To explore the expression 
patterns of the homologous genes between Chinese 
cabbage and Arabidopsis, gene expression data for 
Arabidopsis TCS genes were retrieved from the 
Genevestigator database (data not shown) and then 
we compared these data with our results. Most of the 
homologous genes showed similar expression patterns, 
suggesting the functional conservation of homologs. 
Comparisons of the tissue-specific expression profiles 
between Chinese cabbage and Arabidopsis TCS genes 
might help in determining the unknown functions of 
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Figure 7. Hierarchical clustering and heat map representation for the 
response patterns to drought and high salinity conditions of TCS 
genes in Chinese cabbage roots. The expression levels of genes 
are presented using fold-change values transformed to Log2 
format compared with control. This figure appears in colour in 
the online version of DNA Research. 



the TCS genes in Chinese cabbage. In Arabidopsis, 
AtCKII is mainly expressed in flowers and participated 
in female gametophyte development.^*"^^ Similarly, 
BrHKl , which is a AtCKII homologous gene, showed 
major transcriptions in flowers and pistils in Chinese 
cabbage. Similar expression profiles may suggest 
similar function for Arabidopsis orthologous TCS genes 
in Chinese cabbage. Based on the high sequence identity 
and similar expression patterns, we proposed that 
BrHKl was also related to female gametophyte develop- 
ment in Chinese cabbage. Consistently, with a loss of 
function mutant line, Ito and Kurata^^ revealed that 
OsHKI , which IsaAtCKH homologous gene in rice, is es- 
sential for rice female gametophyte development. This 
result further strengthened our hypothesis and proved 
the conservation function of CKI1 . Furthermore, 
AtAHK5/CKI2 is mainly expressed in roots^^ and regu- 
lated salt sensitivity and resistance against bacterial 
and fungal infection''*'; in addition, AtAHK5/CKI2 
acted to integrate multiplesignalsvia H2O2 homeostasis 
and is independent of ABA signalling in guard cells.^^ 
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Two homologous genes, BrHI<3 and BrHK4, were ana- 
lysed in Chinese cabbage, and these genes also showed 
major transcriptions in roots. Thus, BrHIO and BrHK4 
were probably involved in stress adaptations in 
Chinese cabbage. 

3.6. Expression profiles of TCS genes under drought 
and salt stresses in Chinese cabbage roots 
In addition to cytokinin signalling, accumulating evi- 
dence provesthatTCSisalso involved in plant adaptation 
to various environmental stresses. Relationships 
between cytokinin signalling and stress responses were 
summarized in several reviews.^^'^ °^ Roots were 
generally believed to be the main organs involved in 
drought and high salinity stresses, thus 26 of the 31 
examined TCS genes that showed relatively high expres- 
sion levels in Chinese cabbage roots were selected; in 
addition, qRT-PCR was used to analyse their expression 
profiles under drought and salt stress conditions. As 
shown in Fig. 7 and Supplementary Figs S7 and S8, 
most of the TCS genes, as a whole, were suppressed in 
drought conditions. In different degrees of drought 
symptoms, expressions of BrHKG, BrHK7, BrHI<8, 
BrHP2-3, BrHP4, BrHP7, BrRR1,BrRR3, BrRRZ, BrRR23, 
BrRR27, and BrRR30 were continually suppressed and 
decreased to ~1 5% of the control level. Transcripts of 
BrHK3, BrHIO, BrHPG, and BrRR8 were also primarily 
suppressed, which subsequently increased again after a 
fall in level ?; however, the transcriptions of these genes 
were still below the basal level in level III, which indicates 
that they were potential negative regulators in drought 
stress response. BrHPS, BrHPI , BrRR19, and BrRR31 
were clustered into one group, and these genes showed 
upregulated expressions in levels II or III. BrHPI was con- 
sistently induced; this finding suggests the positive roles 
of BrHPI in drought stress. For the salt treatment, most 
of the TCS genes also showed decreased expression 
levels; however, the responses of these genes were not 
always the same with drought treatment. BrHI<8 and 
BrRR27 seemed to be unaffected with salt treatment; 
however, these genes were continually suppressed in 
drought conditions. The transcript of BrPHPI , which is 
a pseudo-HPgene, was unaffected bydroughttreatment; 
BrPHPI transcript was suppressed in high salinity condi- 
tions. BrRRS was mainly induced with 2-fold upregulated 
expression levels after 48 h salt treatment; however, this 
gene was suppressed in drought treatment. Likewise, 
BrHPS, BrHPI, BrRR19, and BrRRS 1 were clustered to- 
gether with increased expressions,suggestingtheir posi- 
tive functions under high salinity conditions. 6rR/? 7 9 and 
BrRRS 1 transcripts were stable 8 h after salt treatment; 
however, these genes were continually induced after 
8 h and showed 3.4- and 2.0-fold increased expression 
levels at 48 h. Theoretically, almost all the TCS genes in 
Chinese cabbage could be induced by drought and salt 



stresses because of the presence of stress-inducible cis- 
regulatory elements in their promoter regions. 
However, our results conflicted with our hypothesis 
with most of the genes being supressed. Although cis- 
element can be used to predict stress-responsive genes 
via a c/s-ele me nt- based targeted gene-finding ap- 
proach;^^ however, the frequency of a c/s-element se- 
quence in the whole genome is relatively high because 
of theirshort length (5-9 bpcore) and flexibility. In add- 
ition, a number of c/s-element sequences might be syn- 
tactically correct without providing practical regulatory 
function. Moreover, there are also stress- repressive 
genes among the TCS genes and we have limited knowl- 
edge on the stress-repressive c/s-element' °^ In fact, the 
patterns of responses to abiotic stress stimuli of the TCS 
genes were complicated and varied in terms of the 
species, genotypes, and organs or tissues. \nArabidopsis, 
drought significantly induced the expression of a Type- 
AARR gene subset, APRS, ARR7, and ARRIS,^ "'^ whereas 
almost all genes for Type-A RR genes were suppressed 
by drought stress in rice.^^ Moreover, all of the 
members of the TCS family had differential transcript 
abundance (under both non-stress and salinity stress 
conditions) a mongsalt-tolerant (Pokkali) and salt-sensi- 
tive (IR64) rice genotypes.^^ Furthermore, Argueso 
etal.^°^ collected data from BIOArray Resource and illu- 
strated that TCS gene responsiveness to different abiotic 
stresses \nArabidopsis root and shoot tissues are differ- 
ent. Le et a/. ^ ^ exa m i ned a nd com pa red t he TCS gene ex- 
pression profiles in root and shoot tissues of soybean 
plants under dehydration stress; in addition, Le et al. 
found that the majority of soybean TCS genes respond 
to dehydration stress in either root- or shoot-specific 
manner. This is also the case for the expression patterns 
of GmARFs in drought stress.' °^ So c/s-element-based 
predicted genes should be validated experimentally to 
reveal their authentic expression profiles in specific 
tissues or organs upon various abiotic stress conditions. 
One of our interests is that whether genes located in 
the same clade of the plylogenetic tree would show 
similar response patterns to abotic stresses. Thus, we 
carefully examined and compared the expression pat- 
terns of TCS genes in Arabidopsis, rice, and soybean 
with our results. It was amazing that even the homolo- 
gous TCS genes showed varied response patterns, sug- 
gesting the stress-responsive functional divergence of 
the homologous genes. In Arabidopsis, AHK2, AHKS, 
and CRE1 transcripts were all rapidly induced by dehy- 
dration;/\/-//<'2 expression alsoappeared to be influenced 
by NaCi and ABA treatments; and induction of the AHKS 
mRNAwas observed during high salinity and cold stres- 
ses, whereas expressions of CKIl and CKI2 were not 
induced by any treatment.^^ However, in rice, the HK 
genes (HK5 and HKS) that a re homologous to CK recep- 
tor genes (AHI<2,AHI<S, and CREl) \n Arabidopsis were 
induced by drought stress, while HK6 homologous to 
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CKI2 was repressed when exposed to drought stress.^ ^ 
Likewise, the soybean HK genes (GmHKI 0-1 7) that 
are homologous to Arabidopsis CK receptor genes 
(AH l<2, AH l<3, and CRE1) showed more diverse response 
patterns, for example, GmHKI 1 was repressed while 
GmHI<12 was induced under dehydration stress.^^ 
Furthermore, response patterns unlikely correlated 
well with their biological functions under abiotic stress 
conditions. We have mentioned that transcripts of 
AHK2,AHK3, and CRE1 \n Arabidopsis cou\d he induced 
by abiotic stresses;^^ however, loss-of-f unction analysis 
of ahl<2,ahk3, and crel indicated that the stress-respon- 
sive /4H/<'2,/4/-//<'3, and CRE1 acted as negative regulators 
in ABAsignalling; in addition,/4/-//<2 andAHIO negatively 
control osmotic stress responses \n Arabidopsis. Other 
supporting evidence was that AHK2 and AHK3 have 
negative regulatory functions in cold stress signalling, 
via ABA response inhibition; however, the transcripts of 
these genes were not a Itered by cold cond ition.^ ^ The i n- 
formation regarding TCS gene identification would be 
more valuable if the functional predication is verified. 
For the TCS genes in Chinese cabbage, expression profil- 
ing is an initial step; further elucidating the potential 
functionsof these genes involved in abioticstress adapta- 
tions would be an interesting and meaningful study. 



3.7. Effects of exogenous tZ and ABA on TCS gene 
expressions in Chinese cabbage leaves 
The response patterns of TCS genes to exogenous 
cytokinin levels were carefully examined \n Arabidopsis. 
Previous data based on Northern blot, microarray, and 
qRT-PCR analysis indicated that AHK, AHP, and Type-B 
ARR genes were not regulated by exogenous tZ and BA 
except for the upregulated AHK1 and AHK4; Type-A 
ARR genes, which function as the cytokinin primary re- 
sponse genes, are strongly induced upon treatment 
with cytokinins.^°^~' °^ To determine the effects of ex- 
ogenous cytokinin on TCS gene expressions, we investi- 
gated the expression profiles of 2 3 of the 31 examined 
TCS genes with relatively high expression levels in 
Chinese cabbage leaves by qRT-PCR analyses (Fig. 8 
and Supplementary Fig. S9). Spraying is a simple and 
easy method for treatment and also leaves response 
well to exogenous phytohormones.Our results revealed 
that the selected seven Type-A RR genes were transcrip- 
tionally upregulated by tZ.The inductionsof these genes 
occurred in 30 min; the transcripts of such genes 
remained at high levels at 1, 2, and 4 h. Many BrHK, 
BrHP, and Type-B BrRR genes were primarily suppressed; 
afterward, their transcripts increased again and 
remained at the basal levels at 2 and 4 h. Considering 
that Type-A ARR genes negatively regulated HP and 
Type-B ARR gene transcriptions in Arabidopsis,'^^'^ ^ ° we 
supposed that the instantly upregulated Type-A BrRR 
genes suppressed the BrHP and Type-B BrRR gene 
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Figures. Hierarchical clustering and heat map representation for the 
response patterns to exogenous tZ of TCS genes in Chinese cabbage 
leaves. The expression levels of genes are presented using fold- 
change values transformed to Log2 format compared with 
control. This figure appears in colour in the online version of 
DNA Research. 

transcriptions; therefore, these events led to a negative 
feedback loop, thereby suppressing the expressions of 
putative cytokinin receptors {BrHI<6-9) in Chinese 
cabbage; furthermore, autoregulations of these genes 
ensured transcript abundance as time progresses. ABA 
is the key hormone involved in the regulation of 
both stress- and non-stress-related processes.' ^ ^'"^ 
Furthermore, cytokinins are generally postulated to 
function as antagonists of ABA in various growth and 
physiological processes, including environmental stress 
responses. Studies have also suggested the presence of 
intensive interactions and crosstalk between cytokinins 
and ABAaswell as theirsignalling pathways.' °°Thus,we 
also examined the effects of exogenous ABA on TCS gene 
expressions in Chinese cabbage. Figure 9 and 
Supplementary Fig. SI 0 showed that many TCS genes 
were suppressed by ABA treatment; this finding was con- 
sistent with the response patterns of drought and salt 
treatment. BrHK9, BrRR3, and BrRRl 7 were strongly 
and instantly suppressed, and even dropped to un- 
detectable levels after 30 min. BrHP4, BrHP7, BrRR28, 
and BrRR30 transcripts were relatively stable with ABA 
treatment. BrHPl, BrHP6, BrRR 7 , and BrRR32 were clus- 
tered together with increased expression levels upon 
ABA treatment. BrHPl was also induced by drought 
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Figure 9. Hierarchical clustering and heat map representation for the 
response patterns to exogenous ABA of TCS genes in Chinese 
cabbage leaves. The expression levels of genes are presented 
using fold-change values transformed to Log2 format compared 
with the control. This figure appears in colour in the online 
version of DNA Research. 

and salttreatnnent;this result furtherconfirmsthe posi- 
tive regulations of BrHPl in stress environments. 
However,similartothecasewithdroughtandsalttreat- 
ment, promoter regions of Chinese cabbage TCS genes 
were rich in c/s-element involved in ABA response, 
such as MYBR, MYCR, and ABRE etc., many of the TCS 
genes were contrarily suppressed by ABA treatment, 
and these could be mainly attributed to the complex 
interactions of hormone network in addition to the 
own limitations of c/s-element- based predicting ap- 
proach. 



3.8. Conclusions 

In summary, a total of 20 HK(L), 8 HP,and 57 RR pro- 
teins were identified and annotated in the Chinese 
cabbage genome. This study provided the first insight 
into the TCS members of Chinese cabbage. The analysis 
resultsofduplications,evolutionarypatterns,and diver- 
gence of TCS genes gained useful information on the 
evolutionary aspects of Chinese cabbage genome. We 
also focused on the response patterns of the TCS 
genes to drought and high salinity conditions and 
screened numerous candidate stress-responsive genes 
in Chinese cabbage. TCS has an important function in 



signal transduction involved in plant growth, develop- 
ment, and environmental stimuli. Moreover, Chinese 
cabbage is one of the most important vegetables that 
are widely cultivated. Our results contributed relevant 
information to molecular genetic studies, thereby pro- 
viding a better understanding of the biological func- 
tions of the TCS genes in Chinese cabbage. Our study 
may also provide guidance for molecular breeders to 
develop economically important high-yielding and 
high-quality stress-tolerant crops in agriculture. 

Supplementary data: Supplementary data are 
available at www.dnaresearch.oxfordjournals.org. 
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